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Abstract 

The sensitivity of LISA depends on the suppression of several noise sources; dominant 
one is laser frequency noise. It has been shown that the six Doppler data streams obtained 
from three space-crafts can be appropriately time delayed and optimally combined to cancel 
this laser frequency noise. We show that the optimal data combinations when operated in 
a network mode improves the sensitivity over Michelson ranging from 40% to 100%. In this 
article, we summarize these results. We further show that the residual laser noise in the 
optimal data combination due to typical arm-length inaccuracy of 200 m is much below the 
level of optical path and the proof mass noises. 



1 Introduction 

The future space-based gravitational wave (GW) mission^- the Laser Interferometrie Space 
Antenna (LISA) - consists of three identical space-crafts forming an equilateral triangle of 
side 5 x 10 6 km following heliocentric orbits trailing the Earth by 20°. LISA is thus a giant 
interferometrie configuration with three arms which will give independent information on 
GW polarizations and will detect GW in the low frequency range of 0.1 mHz to 0.1 Hz. 
Due to the long arm-lengths of the antenna, it is not feasible to bounce the laser beams. 
A special Doppler tracking scheme is used to track the space-crafts with laser beams. This 
exchange of laser beams between the three space-crafts result in six Doppler data streams. 
The LISA sensitivity is limited by many noise sources; the dominant one is the laser phase 
noise; noise due to phase fluctuations of the master laser. The current stabilization schemes 
estimate this noise to about AvjvQ ~ 10~ 13 /VHz, where Vq is the frequency of the laser 
and Av the fluctuation in frequency If the laser frequency noise can be suppressed then 
the noise floor is determined by the optical-path noise which fakes the fluctuations in the 
lengths of optical paths and the residual acceleration of proof masses resulting from imperfect 
shielding of the drag-free system. Thus, canceling the laser frequency noise is vital for LISA 



Figure 1: The LISA constellation 



to reach the requisite sensitivity of h ~ 10 21 or 10 22 . Since it is impossible to maintain 
equal distances between space-crafts, cancellation of laser frequency noise is a non-tr ivial 
problem. Several schemes have been proposed to combat this noise. In these schemes^El 
the data streams are combined with appropriate time delays in order to cancel the laser 
frequency noise. In our earlier work, it was rigorously shown that all laser-noise free data 
combinations form an algebraic module over a polynomial ring over time delay operators^. 
Furthermore, recently^! we have found data combinations which are eigen-data combinations 
of the noise covariance matrix, which give optimal sensitivity, when averaged over all the 
source directions and polarization angles. We show that signal-to-noise ratio (SNR) for any 
data combination in the laser noise-free cancellation module lies between th e SN R's of these 
eigen-data combinations. In this article, we summarize our previous resultsEEl. Besides, we 
estimate the residual laser noise due to inaccuracies in the arm-lengths. 



2 Time Delay Interferometry and Laser Noise Cancellation Module 

In LISA constellation, the six Doppler data streams are labeled as U % and V 1 , i = 1,2,3; if 
the space-crafts are labeled clockwise as shown in Fig. ^ The Doppler data stream U 1 is 
obtained by letting the laser beam from space-craft 3 to travel towards space-craft 1 along 
the arm of length L 2 in the direction — A2, and is beaten with the on-board laser beam at 
space-craft 1. Similarly, — V 1 represents the Doppler data obtained by beating laser beam 
traveling from space-craft 2 to space-craft 1 along the arm of length L3 in the direction of A3 
with the on-board laser at space-craft 1. The remaining 4 beams are described by cyclically 
permuting the indices. These beams contain the laser frequency noise, other noises such as 
optical path, acceleration etc. and also the GW signal. Using time delay interferometry, 
these 6 Doppler data streams can be appropriately combined by delaying them with time 
delay operators Eix{t) = x(t — L{),i = 1,2,3 such that the resultant data streams is laser 
noise-free. Any laser noise-free data combination A is represented by 

A = Pt V* + q i U i , (1) 

where pi, qi are polynomials of time-delay operators Ei. For the LISA configuration Li ~ 16.7 
seconds, corresponding to an arm-length of 5 million km, one such known laser noise-free 
data combination, expressed in terms of 6-tuples of polynomials (pj, qi) is 

X = (1 - El 0, E 2 (El - 1), 1 - El,E 3 (E 2 2 - 1), 0) (2) 

commonly referred to as Michelson combination in the literature^ in which one arm of 
LISA is not used. In general, we have shown^ that all the data combinations which cancel 



the laser frequency noise and the optical bench motion noise form an algebraic module of 
syzygies over a ring of polynomials in time delay operators Ei . This formalism generates the 
noise cancellation module from the generators; linear combinations of the generators with 
polynomial coefficients in the ring generates a module. One such set of generators (whic h is 
convenient for SNR optimization purpose) is a,/3, 7 and £ (notation followed fronJ * ' j) . 

a = (1, E 3 , EiE 3 , 1, EiE 2 , E 2 ) , 

f3 = (EiE 2 , 1, Ei, E 3 , 1, E 2 E 3 ) , 

7 = (E 2 ,E 2 E 3 ,l,E 3 E 1 ,Ei,l), 

C = (Ei, E 2 , E 3 , Ei, E 2 , E 3 ) . (3) 

We note that a, (3, 7 are cyclic permutations of each other. The combination £ is the 
symmetric Sagnac combination which is insensitive to GW at low frequency due to its high 
symmetry . 



3 LISA Sensitivity Optimization 

To generate any laser noise-free data combination, in general, 4 generators are necessary 
but if the source is monochromatic, the fourth generator £ can be effectively eliminated by 
expressing in terms of (a, (3, 7) as follows 

(1 - EiE 2 E 3 )( = {Ei - E 2 E 3 )a + (E 2 - EiE 3 )(3 + [E 3 ~ EiE 2 )~f . (4) 

except at certain frequencies which are solutions of e l ( L i+ L 2+L 3 )fi _ 1 As the maximization 
is possible arbitrarily close to the singular frequencies, the singularities do not seem to be 
important. 

Since the difference in arm-lengths of LISA is smaller than the GW wavelength^, for 
computing the response, all the arms can be taken to be equal. This simplifies further analy- 
sis. We then show that the set of 3 optimal data combinations having noises uncorrelated to 
each other can be obtained by linearly combining a,f3,"f. This new set acts as optimal when 
one averages over all directions and polarizations of the binary system. Below, we briefly 
summarize this optimization. 

(a) Noise covariance matrix: We define noise vectors in the Fourier 

1,2,3 for each of the generators X^- 1 ' = {a, /3,7}, respectively, over the 12 dimensional 

complex space C 12 , 



where S pi (/) = 2.5 x lfT 48 [f/1 Hz]~ 2 Hz~ l and S opt {f) = 1.8 x 1CT 37 [f/l Hz} 2 Hz' 1 arc 
power spectral densities (psd) of the proof mass residual motion and the optical path noise 
respectively^. The polynomials (pP , qP^) corresponding to the generators are given in 
the equation ©. The rP polynomials are defined through rP = —(pP+E 3 qP ) = — (qP + 
E 2 pP) plus cyclic permutations for r*p and rP . For a given data combination, the norm 
of the noise vector represents its noise psd. The noise covariance matrix Npi — ■ Nfjs 
defined for above generating set X^' takes a symmetric form as 

= {n d for I = J and n D for I ^ J}. (6) 

(b) Signal covariance matrix: The response of a GW signal for a given laser noise-free 

m 

data combination A is expressed in the Fourier domain 1 ^ as, 

3 

h (A \n) = [Pi A) (Fv i ., + h+ + Fy^hy) + qP ] (F Uh+ h + + F v . x h x )] (fi) . (7) 
i=i 
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Figure 2: Log Log plot of sensitivity S, curve as function of / after averaging over polarization and source 
directions for a observation period of one year and SNR =5. 



Here, F Vi . +/x a,mi r, 



and Ftl +/x are the antenna pattern functions. For a binary source which may 



be adiabatically changing in frequency, the two GW amplitudes at frequency f2 are given by, 



1 + cos^ e 
2 

1 + cos 2 e 



cos 2-0 — i cos e sin2ij} 
sin2ip — i cos e cos 2ip 



h+(Sl) = A 
h x (Q) = A 



Here, the polarization angles e and ip describe the orientation of the source and enter into 
the expressions for the polarization amplitudes. The direction of the source on the celestial 
sphere is given by the angles 9 and <j>. Similar to noise covariance matrix, the signal covariance 
matrix averaged over all directions and polarizations is defined as, 



(8) 



where ( ) e ^e4> represents the average over the polarizations and directions. We note that 
T6 IJ } takes the same symmetric structure as the noise covariance matrix given m 

equation JfJJl. This is due to inbuilt cyclic symmetry in this generating set Thus, the 

diagonal component of liS 1 ^ is hd and the off-diagonal is h . 

Due to above symmetric structure, both liS 1 ^ and M^ 1 ^ can be diagonalized simulta- 
neously. The common eigen-observables thus obtained are given by 



y(i) 



(a + /? - 2 7 ) , = _L _ a ) , yO) 



v/3 



(a + (3 + 7) 



(9) 



V6 y " ' V2 

The above eigen-observables have following properties: 

• Data combinations and Y^> have same SNR given by \l (hd — h )/ (rid — n ) ■ 
The SNR of Y^ is yj (h d + 2h ) / (n d + 2n ) . At any given frequency, if SNR of rW 
is greater than SNR of Y^ 3 \ then Y^ gives maximimum SNR whereas Y^ gives 
minimum SNR or vice-a-versa. Thus, at any given frequency, either Y"W or Y^ acts 
as optimal data combination amongst all the data combinations in the laser noise-free 
module. 

• Both and Y^ perform comparable to the Michelson combination X at low fre- 
quency. Whereas Y^ is proportional to the symmetric Sagnac and hence insensitive 
to GW at low frequencies (~ below 3 mHz). 

We note that averaging over directions and polarizations^ results in signal covariance matrix 
of rank 3. We note that if we average over the polarizations only and obtain the optimal 
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Figure 3: Plots showing the relative improvements (ratios) of SNRs for the three cases: (i) Network SNR over the 
Michelson data combination (solid line), (ii) Network SNR over the maximum of Max [yW-yC 3 )] (dotted line), 
(iii) Max [yM-yW] over the Michelson (dashed line). Here Max [Y w 'Y i3) ] is the maximum of the SNR of Y w 

and y (3) over the bandwidth of LISA. 



data combination for a particular direction, the signal covariance matrix is of rank 2. Such 
an optimal combination is important while optimally tracking the source in LISA framed 
We may further note that in the formalism developed by Prince et aJ^j the optimization is 
performed without averaging over the source directions and polarizations, which results in 
signal covariance matrix of rank 1. 

We have shown that either or maximize the LISA sensitivity on an average 

sense and they are orthogonal i.e. they are independent random variables. The sensitivity 
of LISA can be further improved as each of these generators can be realized as independent 
gravitational wave detectors. We assume that the Y^ follow the Gaussian noise distribution 
and thus quadratically can combine the SNR's of these eigen-observables to form a network- 
observable^. The network SNR is given by 

3 

SNRl etwork = SNR h = 1SNRl (1) + SNR Y{3) . (10) 

7=1 

The corresponding sensitivities are shown in the Fig. [3 In Fig. |3| we have plotted the 
relative improvements in the network SNR with respect to the Michelson combination and 
the maximum of Y^ and Y ^ . At low frequencies / < 15 mHz, the improvement of the 
network SNR over the maximum of Y^ 1 ' is slightly greater than v2< This is because at low 
frequencies the data combination Y^ is not very sensitive in comparison with yW. The 
best improvement of factor v3 in the relative SNR is achieved at frequencies where all the 
data combinations are equally sensitive, that is, when SNRy(i) = SNRy(3). 



4 Residual laser noise 



As we know, it is difficult to maintain constant distance between the three space-crafts. In 
general, the three arm-lengths will be different. The actual length L[ can be estimated upto 
a certain accuracy. Let the estimated arm-length be Li with an error ALi such that the 
actual length is L\ = L% + ALi. Because of this unknown inaccuracy, the laser noise is not 
completely canceled in a given laser noise-free data stream. If we demand that this residual 
laser phase noise level should be below the combined noise from proof mass S p f and the 
optical path noise S op t then this puts an upper limit on the laser stabilization requirement 
Av and is given by 



Av = 



Spf + S, 



nAL [J2\ Pz 



opt 



1/2 



(11) 




1 o' 



10 4 10 3 

Frequency in Hz 

Figure 4: Laser frequency stabilization Av in Hz/-\/Hz as function of frequency for Y' 3 ' and Michelson X 

combinations, for AL — 200 m. 



For Nd-YAG laser, v = 3 x 10 14 Hz and the inaccuracy in length is assumed to be AL ~ 200 
m, the upper limit on the laser stabilization requirement is plotted for data combinations Y/n 
in Fig. For this assumed inaccuracy in arms, the combination demands that the laser 
frequency stabilization be at least as good as Av ~ 25 Hz/Viz" at 1 mHz. While for Y {1 \ 
the requirement is much less stringent on frequency stabilization. The laser stabilization 
requirement scales linearly with the assumed arm-length inaccuracy. 

5 Conclusion 

For any frequency in LISA band, we show that the optimal laser noise-free data combina- 
tions (when averaged over all directions and polarizations of GW source) are nothing but 
the eigen-data combinations of noise-covariance matrix. Since these combinations have un- 
correlated noise, their SNR's can be combined quadratically to improve LISA sensitivity. 
The improvement varies from 40% to 100% with respect to Michelson data combination. 
We further show that with our demand that the residual laser noise should be less than 
the proof mass noise and optical path noise, the laser frequency stabilization requirement 
varies inversely proportional to the arm-length inaccuracy. The stringent demand on the 
laser stabilization requirement is Av ~ 25 Hz/VHz at 1 mHz for inaccuracies in arm- length 
of 200 m. 
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